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Abstract 

We present results on the system size dependence of high transverse momentum 
di-hadron correlations at ^/snn = 200 GeV as measured by STAR at RHIC. 
Measurements in d+Au, Cu+Cu and Au+Au collisions reveal similar jet-like 
correlation yields at small angular separation {Acj) ^ 0, Ar] ^ 0) for all systems 
and centralities. Previous measurements have shown that the away-side yield is 
suppressed in heavy-ion collisions. We present measurements of the away-side 
suppression as a function of transverse momentum and centrality in Cu-|-Cu 
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and Au+Au collisions. The suppression is found to be similar in Cu+Cu and 
Au+Au collisions at a similar number of participants. The results are com- 
pared to theoretical calculations based on the parton quenching model and the 
modified fragmentation model. The observed diff'erences between data and the- 
ory indicate that the correlated yields presented here will provide important 
constraints on medium density profile and energy loss model parameters. 

Key words: parton energy loss, jet quenching, di-hadron fragmentation 
function, relativistic heavy-ion collisions 
PACS: 25.75.-q, 25.75.Gz 



One of the important early results from the experiments at RHIC is the 
observation of jet quenching in heavy-ion collisions. The suppression of high 
transverse momentum (px) particle production in inclusive hadron spectra 
and jet-like structures in di-hadron correlation measurements Q indicates that 
partons originating from hard scatterings in the initial stages of the collisions 
interact strongly with the created medium. 

Previous studies Q investigated azimuthal correlations of high-p^ hadrons 
and showed that the suppression of the correlated away-side yield increases with 
centrality in Au-f-Au collisions. Various theoretical calculations [^, @, 0, 0] of 
partonic energy loss in the medium have been used to derive medium prop- 
erties like the transport coefficient q. The energy loss of individual partons 
is also expected to depend on the path length through the medium in a way 
that is characteristic of the energy loss mechanism. For radiative energy loss, 
which is thought to be dominant for light quarks, the energy loss is expected 
to depend on L^, the square of the traversed path length, due to coherence 
effects 

[i, Eini- For elastic energy loss, on the other hand, a linear depen- 
dence on L is expected. Prior results from RHIC have not established in detail 
the energy loss mechanism. Combined measurements of single-hadron and di- 
hadron suppression are sensitive to the path length dependence and can help 
determine which process dominates 12| . In addition, different implementations 
of the energy loss calculation use different path-length distributions and density 
profiles. The system-size dependence of away-side suppression is sensitive to 
these modeling parameters and will provide further constraints [l3| . 

We present a systematic study of the near- (Acj) ~ 0) and away-side {\A(j)\ ~ 
tt) di-hadron correlated yields as a function of the number of participant nucleons 
impart)- Data for three systems with diflFerent geometries (d+Au, Cu-|-Cu and 
Au-I-Au) at y^sWN = 200 GeV were collected by the STAR experiment at RHIC. 
A study of the hadron-triggered fragmentation functions in the three systems 
is also presented. Results from d-|-Au collisions are used as a reference without 
a hot medium. The d-|-Au data sample is preferred over the p-l-p data sample 
because it has significantly larger statistics. Earlier comparisons between p-l-p 
and d-|-Au collisions have established that jet suppression is a final state effect 
and is not present in d-|-Au collisions [14, 15, H, 17 1. 

This analysis is based on four data sets and includes 11.7 million minimum- 
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Figure 1: Di-hadron correlations in central (0-12%) Au+Au collisions: (a) Acj) 
correlations - small Arj {\Ari\ < 0.7) (black circles) and large At; (0.7 < \Ar]\ < 
1.7) scaled to match the small A77 result at large Acf) (red triangles), (b) Atp 
subtracted distribution, (c) A77 subtracted distributions; 4 < p^'^ < 6 GeV/c, 
pf""'' >3 GeV/c. 
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bias d+Au events, 43.8 million minimum-bias Cu+Cu events, 25 million minimum- 
bias Au-|-Au events and 19 million Au-|-Au events collected using a central trig- 
ger. The central trigger uses the coincidence of two Zero Degree Calorimeters 
(ZDCs) and a multiplicity threshold in the Central Trigger Barrel [31 which 
selects the most central — 12% of total geometric cross-section. In order to 
minimize the influence of the Time Projection Chamber (TPC) acceptance [l9j . 
only events with a reconstructed collision vertex with \zvertex\ < 30 cm are in- 
cluded in the analysis, where z^ertex represents the distance along the beam-line 
from the center of the detector. 

The di-hadron correlations are formed using charged particles reconstructed 
in the TPC, within a pseudorapidity range of —1 < 77 < 1. High pT trig- 
ger particles are selected and the A77 x A</) distribution of associated particles 
^p^assoc ^ p^'^) is constructed. An 77, pT and centrality dependent recon- 
struction efficiency correction is applied to obtain the associated particle yields. 
It is not necessary to apply the efficiency correction to trigger particles when 
calculating the correlated yields because the final result is normalised per trig- 
ger particle. The track reconstruction efficiency depends on the track density 
within the TPC and ranges from 89% (peripheral collisions) to 77% (central). 
The systematic uncertainty on the efficiency correction is estimated to be 5% 
and is strongly correlated across centralities and pr bins for each data set, but 
not between data sets, particles is not necessary since 

Pair acceptance corrections in A(f> and A77 are computed using a mixed event 
technique. These corrections reflect the conditional probability of reconstructing 
two tracks with a specified relative kinematics. The dominant feature in the Acj} 
pair acceptance correction are the small gaps between the sectors of the TPC. 
The A77 pair acceptance correction is of triangular shape, with a maximum of 1 
at mid-rapidity and minimum of at the limit of the pair acceptance Ar; = ±2. 
The A77 acceptance correction is not applied to the away-side yields. 

Earlier results from STAR [2^, [2l[ have shown that there is a finite associ- 
ated yield on the near-side {A(j) ~ 0) with large pseudorapidity separation A77 
(the "ridge"). Since the ridge properties are similar to those of the medium, 
it is appropriate to subtract this contribution in the present analysis. In order 
to extract the jet contribution to the near-side yield, the azimuthal correlation 
distribution for large At] separation (0.7 < IA77I < 1.7) is subtracted from the 
distribution for small At] (lAjyj < 0.7). To account for the different At] window 
widths, the former distribution is scaled so that the two distributions match 
in the away-side region. This subtraction removes the A77 -independent ridge 
contribution and the contributions from elliptic flow W2, which is also largely 
independent of ry in the range considered [l^l . Figure [Ta| shows central Au-|- Au 
distributions in the large (black) and small (red) A77 regions, for trigger parti- 
cles with transverse momentum 4 GeV/c < pt*"^ < 6 GeV/c, and associated 
particles with prp"-^^°''- in the range 3 GcV/c < pj'°-'^^°^ < prp^"^ . The signal 
distribution after subtraction is shown in Fig. Ilbl 

An alternative way to extract the near-side associated yield is to use the 
Ary-distribution, which is obtained by projecting the Ar] x Acf) correlations in 
the |A^| < 0.78 region onto the Ar] axis. In this projection, the ridge yield 
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Figure 2: Npart dependence of the the near-side associated-particle yield for 

two trigger pT ranges: (a) 4 GeV/c < p^*^ < 6 GeV/c, (b) 6 GeV/c < 

pp^ < 10 GeV/c. For both panels 3 GcV/c < p^''°'= < pp^ . The hollow 
symbols arc horizontally offset for clarity. 



and elliptic flow constitute a flat background which is determined by averaging 
the yield at large lAryj > 0.7 and subtracted. Figure [Tel shows a background 
subtracted Ar? projection with the same trigger conditions as Fig.lTaland Fig.lTbl 
The near-side associated-particle yield, defined as 

M.7 1-0.78 ,2jsj 

J -0.7 7-0.78 Ntrig d{Ari)d{A(j>) 

is presented as a function of number of participant nuclcons (Npart) in Fig. [2l 
The two methods produce results that arc consistent with each other. The 
Cu-|-Cu and AuH-Au near-side associated yields arc consistent within errors 
for similar Npart- The near-side yields in heavy- ion collisions show no centrality 
dependence and within errors agree with those in d-|-Au, as seen also in previous 
studies 0]. The observed independence of the near-side associated yields on 
centrality indicates that in this pr-range fragmentation is largely unmodified 
by the presence of the medium. Note that this does not necessarily imply 
that those partons do not lose energy, but rather that they fragment outside 
the medium after energy loss. In that case, the energy loss would reduce the 
number of trigger hadrons at a given pT, but not change the associated particle 
distribution at intermediate to high pT- The enhancement of associated particles 
at low pt that has been reported earlier j20| could then be due to fragments of 
radiated gluons. 

The choice of high-pT trigger particles leads to a surface bias in the dis- 
tribution of hard scattering points The away-side partons have longer 
path lengths through the medium and therefore will suffer higher energy losses 
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Figure 3: Acj) distribution in (0-12%) central Au+Au collisions used to extract 
the away-side yield, 4 GeV/c < pp^ < 6 GeV/c, 3 GeV/c < p^f"'"^ < pp^, 
\Ar]\ < 1.7. The triangular pair acceptance correction in Ar; is not applied. The 
green line represents the elliptic flow modulated background using the values of 
calculated using the reaction plane method, the blue line uses the V2 obtained 
using the 4-particle cumulant method. The red line uses the average value of 

V2- 



that lead to away-side yield suppression. The study of the away-side yield 
suppression provides an important tool for determining the energy loss de- 
pendence on path length. The away-side associated-particle yield is measured 
by integrating the associated hadrons in the region \A(f) — 7r| < 1.3, covering 
the azimuthal range of the away-side jet. A background subtraction is ap- 
plied to remove the background which is correlated with the trigger particles 
through elliptic flow V2- The elliptic flow modulated background is described 
by dN/d{A(t)) = s(^l + 2(4"^u5''"°'=)cos(2A0)), and is illustrated in Fig. [3] 
for central collisions. The background is subtracted using the assumption that 
there is no jet contribution at the minimum of the distribution (23l | — in this 
case at |A0| ^ 1. 

The amplitude of the background modulation is given by {v^^ v!^"'^"'^) ~ 
(v^"^) (yS;''^"'^) which is measured in STAR using a number of different meth- 
ods [2J]. For the Au+Au collisions, the nominal value of V2 for the background 
subtraction was the average between the four-particle cumulant and the reac- 
tion plane measurements of V2. In the Cu-|-Cu case, the nominal value is the 
average between the V2 results obtained using two methods. The first method is 
the reaction plane method using tracks in the Forward Time Projection Cham- 
ber [2^. The second method uses tracks in the TPC but subtracts the azimuthal 
correlations in p+p collisions to remove non-flow correlations. The systematic 
uncertainty associated with the background removal is estimated in both cases 
as the difference between the results given by each method and the nominal 
value. 

The background subtracted away-side yields are used to compute the sup- 
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Figure 4: Npart dependence of the away-side associated-particle yield for 
two trigger px ranges: (a) 4 GeV/c < p*P^ < 6 GeV/c (b) 6 GeV/c < 
pp^ < 10 GeV/c. For both panels 3 GeV/c < p^f""^ < pp^ . The error bars 
represent statistical errors and the boxes represent the point-to-point systematic 
errors. The gray band represents the correlated error due to the statistical error 
in the d+Au data. The lines represent calculations in PQM and MFM models. 
The values of q are expressed in GeV^/fm. 



pression factor Iaa = ^aa^^ l^dAu^ ^ where yXA{dAu) away-side di-hadron 

correlation strength in heavy-ion and d-|-Au collisions, respectively. Figure H] 
shows the results for Iaa as a function of number of participants for Cu-|-Cu 
and Au-|-Au collisions. The away-side yield suppression increases with Npart, 
as expected. The Cu-I-Cu results show a similar suppression (Iaa) at the same 
number of participants as the Au+Au results, despite possible differences in 
density and path length distributions. 

Figure [5] also shows two model calculations implementing the same kine- 
matic cuts as our analysis^ One calculation, the Parton Quenching Model 
(PQM) 0, E^, uses the Salgado- Wiedemann quenching weights (l^l with a 
Glauber-overlap geometry in which the local density scales with the local den- 
sity of binary collisions pcoii ■ The other model uses a next-to-leading order QCD 
calculation with modified fragmentation functions from a higher-twist formal- 
ism [2^ and a hard-sphere geometry where the density scales with the local 
participant density ppart [ll- We refer to this model as the Modified Fragmen- 
tation Model (MFM). The MFM authors used previous data on the suppression 
of high-pT away-side yields in central Au-|-Au collisions [3| to tunc their model. 
The PQM authors present 3 calculations, based on 3 values of q in central 



^The model calculations use p+p as the reference, which is expected to be equivalent to 
the d+Au measurement used in the data. 
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Figure 5: Away-side associated particle distribution and Iaa for 6 < 
p^'^ <10 GeV/c. The error bars represent statistical errors and the boxes 
represent the total systematic errors. The lines represent calculations in MFM 
model. 



collisions, indicated by different line styles in the figure. 

For the lower trigger selection, 4 GcV/c < p^'^ < 6 GcV/c, the Modified 
Fragmentation Model predicts a smaller suppression than observed in the data, 
whereas PQM cannot explain Cu+Cu or Au+Au results in a consistent fashion. 
The disagreement between the models and the data suggests that the effect 
of kinematic limits (energy loss cannot be larger than the jet energy) and non- 
perturbative effects, which are not explicitly treated in the model, are significant 
in this pT-range. For the higher trigger pT range, 6 GeV/c < p^^^ < 10 GeV/c, 
a better agreement between the data and MFM is observed. There is an obvious 
difference between the system size dependence in the two models. While MFM 
obtains Iaa values that are independent on the system at a certain Npart, 
PQM shows a clear difference between the two systems for similar Npart, when 
using a common scaling of the medium density (represented by line styles in 
the figure). Further model studies are needed to clarify whether the different 
scaling behavior in MFM and PQM is mainly a result of the different quenching 
formalisms or rather due to differences between the medium density models. 

In Figs. [2] and |4] we have presented results for a single selection of asso- 
ciated liadrons, p!^**""^ > 3 GeV/c. A more differential measurement is pre- 
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sented in Fig. which shows the away-side associated yield as a function of 
zt = p^^°'^ /p*P^ ■ The lower panel of Fig. [5] shows the ZT-dependence of Iaa- 
The away-side suppression is approximately independent of zt in the measured 
range, indicating that the momentum distribution of fragments along the jet axis 
is not modified by energy loss. A possible explanation of the ZT-independent 
Iaa is that energy loss is large enough that partons which lose energy have 
such a soft fragment distribution that they do not contribute significantly to 
the away-side yield. The remaining away-side yield would then be dominantly 
from the fraction of partons that lost little or no energy due to a short path 
length (surface bias, tangential jets) or energy loss fluctuations. Also shown in 
Fig. O are calculations in the Modified Fragmentation Model which agree 
with the results within the present statistical uncertainties. 

In summary, we have presented a systematic study of di-hadron correlations 
of particles associated with high transverse momentum trigger hadrons. We 
have studied the jet-like correlations on the near-side (A^ ~ 0) and away-side 
(A(/) '-^ tt) for d-l-Au, Cu-I-Cu and Au-|-Au collisions at ^sj^n = 200 GcV/c. 
Near-side associated yields are equal within the experimental uncertainty for 
all the systems studied and independent of the number of participant nuclcons 
(Npart)- Away-side associated yields are suppressed in heavy- ion collisions with 
respect to the d-|-Au reference. The suppression increases with increasing Npart 
and shows no significant dependence on the collision system for a given Npart- 
The Parton Quenching Model does not describe the similarity of the away- 

side yields in the two collision systems at a given Npart, while the Modified 
Fragmentation Model 0, [2^ describes this relatively well for the higher pT 
triggers. Further comparison of these measurements to models may allow the 
extraction of the path length dependence of energy loss and whether elastic or 



radiative energy loss is dominant [12 1 
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